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Within the framework of boundary conformal field theory, we evaluate the conductance of stable fixed points 
of junctions of two and three quantum wires with different Luttinger parameters. For two wires, the physical 
properties are governed by a single effective Luttinger parameter for each of the charge and spin sectors. We 
present numerical density-matrix-renormalization-group calculations of the conductance of a junction of two 
chains of interacting spinless fermions with different interaction strengths, obtained using a recently developed 
method [Phys. Rev. Lett. 105, 226803 (2010)]. The numerical results show very good agreement with the 
analytical predictions. For three spinless wires (i.e., a Y junction) we analytically determine the full phase 
diagram, and compute all fixed-point conductances as a function of the three Luttinger parameters. 

PACS numbers: 73.63.Nm, 71.10.Pm, 73.23.b 



I. INTRODUCTION 

Conducting quantum wires, at low energies, generically 
form a Tomonaga-Luttinger liquid (TLL), characterized by a 
Luttinger parameter g, which encodes the effects of electron- 
electron interactions. A- Due to the prominent role of interac- 
tions in low dimensions, the nature of these one-dimensional 
electronic systems is dramatically different from their higher- 
dimensional counterparts described by Landau's Fermi-liquid 
theory A The TLL state of matter in one-dimensional quantum 
wires has been realized in numerous experiments over the last 
few years«^ 

Transport properties of such quantum wires are of consid- 
erable interest: From a fundamental point of view, a large 
number of interesting phenomena have been predicted and ob- 
served. For instance, at low temperature and low bias voltage, 
a TLL with repulsive interactions (g < 1) is totally discon- 
nected in the presence of an impurity, while one with attractive 
interactions {g > 1) conducts as in the absence of the impu- 
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From a practical viewpoint, junctions of TLL wires 



serve as important building blocks of quantum circuits.^ 2 ! 
and are thus of technological significance. Junctions of three 
quantum wires, known as Y junctions, also have highly non- 
trivial transport properties ^122 Due to their rich transport be- 
havior, junctions of quantum wires and their networks have 
thus attracted much attention. 22r2L 

Most of the previous works on the transport properties of 
junctions of TLL wires focus on wires with the same Lut- 
tinger parameter. However, experimentally, there is no reason 
for all the TLLs emanating from a junction to be identical. 
Moreover, a single TLL can have inhomogeneities; for exam- 
ple, a contact between an interacting TLL and a Fermi-liquid 
lead, a key ingredient of most transport measurements, is of- 
ten studied as an inhomogeneous TLL wire smoothly interpo- 
lating between interacting (TLL) and noninteracting (Fermi- 
liquid) regions or as a two-wire junction with the Luttinger pa- 
rameter abruptly changing at the junctionals A junction of 
three quantum wires with different Luttinger parameters has 
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FIG. 1: Junctions of two and three TLLs with different Luttinger 
parameters. 



been studied in the weak coupling regime^^r— The exper- 
imental importance of junctions of TLL wires with generally 
unequal Luttinger parameters motivates an in-depth study of 
their properties, which is the main objective of the present pa- 
per. 

Here, we focus on junctions of two and three nonchiral 
Luttinger liquids schematically depicted in Fig. Q] For two 
wires, it is known that the transport properties of the junc- 
tion are fully controlled by one effective Luttinger parameter 
g e = 2/igT 1 + g, 1 ) as found in Ref. [39l In the context of 
fractional Hall edge states^, similar results have been found 
for tunneling between two chiral-TLL edge states. 43 . For two 
nonchiral wires, one can reach the same conclusion through 
an almost identical argument. In this work, however, we ob- 
tain this result within the framework of boundary conformal 
field theory (BCFT), using the delayed evaluation of bound- 
ary conditions (DEBC) method^— which, as we will see, 
has the advantage that it can be readily generalized to junc- 
tions of more quantum wires. 

Such generalization to a junction of three nonidentical 
quantum wires is a key result of this paper. We find the sta- 
bility regions of the previously identified (in Ref. [27| for three 
equal Luttinger parameters) fixed points of such a Y junction 
in the (gi,g2, gi) space, and compute their corresponding con- 
ductances as a function of these three Luttinger parameters. 
Moreover, we obtain new asymmetric fixed points, which are 
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only stable for nonidentical TLLs, thereby providing a more 
complete classification of the conformally invariant BCs for 
three TLLs. Such asymmetric fixed points have been identi- 
fied using perturbative renormalization-group analysis in the 
weak coupling regime i 21,57 ' 58 

Another important result of this paper is a direct numerical 
verification, through DMRG computations, 65 of the analytical 
predictions for the conductance of a junction of two noniden- 
tical wires. Using a recently developed method^ 6 ^ which al- 
lows us to extract the conductance from a ground-state calcu- 
lation in a finite system, we compute, in a microscopic lattice 
model, the conductance of a junction of two chains of interact- 
ing spinless electrons with different interaction strengths. Our 
numerical results show excellent agreement with the DEBC 
predictions. 

The outline of this paper is as follows. In Sec. [II] we set up 
the notation and present the model in the bosonization frame- 
work. In Sec. UnJ we present the DEBC analysis of a junc- 
tion of two spinless wires and show that the scaling behav- 
ior of such junctions is governed by a single effective Lut- 
tinger parameter g e = 2/(g~ l + g^ )■ For wires with spin- 
1 /2 electrons, we obtains two effective Luttinger parameters 
g°f = 2/(l/gj' s + l/g")' corresponding to charge and spin 
sectors. Section HVl contains the numerical DMRG calcula- 
tions of the conductance of a junction of two spinless chains 
with different interaction strengths. In Sec.|V] we summarize 
our results on the Y junction with the detailed DEBC analysis 
presented in Appendix lAl In Sec. lVIl we present the analysis 
of the conductance renormalization when the wires are con- 
tacted to Fermi-liquid leads. Finally, we conclude in Sec. lVIII 



II. GENERAL SETUP 



fields, can be written in terms of the above bosons through: 



V2^ 



1> ifa>,-±8,0/ a/2 



V2^ 



(2.4) 



where r\-, are anticommuting Klein factors, which ensure the 
correct fermionic statistics. The Klein factors play no role 
in our analysis and are hence neglected throughout the paper. 
It is convenient to define complex variables z = t + ix and 
Z = r - ix such that the left- and right-moving bosons, corre- 
sponding to current flowing toward and away from the junc- 
tion, are, respectively, functions of z and z only. The chiral 
current operators can then be written as 



J? 



-d-A\ 



V2tt 



V2tt 



-d,6i, 



(2.5) 



where d z = (d T -id x )/2 and 3; = (d T +id x )/2, The total current 
is proportional to the difference between the right and the left 
currents: /, = Vi(jf - Jf). 

To analyze the junctions of quantum wires with unequal 
Luttinger parameters, it is convenient to introduce rescaled 
bosonic fields 



Oi = Oil ^[g~c, (fit = y/gi<Pi, 



(2.6) 



which effectively have a noninteracting action [i.e., g = 1 in 
Eq. d2.1| )l. Note that these rescaled fields satisfy the origi- 
nal commutation relations. Similarly, we define the following 
rescaled left- and right-moving bosonic fields: 



9i = ft + 0f ; 



(2.7) 



in terms of which the left- and right-moving fermions become 



In this section, we present the model in the bosonization 
framework and set up the notation. In the low-energy limit, the 
wires are described by TLLs with the Euclidean action:^^ 



S — V* -^r- I dr I dx 
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where g j and v, are, respectively, the Luttinger parameter and 
the plasmon velocity of wire i. Different wires can have dif- 
ferent electron-electron interactions and consequently differ- 
ent Luttinger parameters g : . The boson fields, <pi and 6>,, have 
the following equal-time commutation relation, 



[<fii(x), Oj(x')] = in6jjSgn(x' - x), 



(2.2) 



so the conjugate momenta of fields are given by lie, = 
(d x (fj)/27T. Let us also define left- and right-moving boson 
fields as 



(2.3) 



The fermions if/i(x) = e' kFX ipf(x) + e ikFX i/ff(x), with ifrf(x) 
and (AfW the linearized left- and right-moving fermionic 



<Af - R = exp 



V2 



(2.8) 



In the absence of a junction (boundary), the correlation func- 
tions of the rescaled fields are given by: 



(<Pi(z> Z)<Pj(w, w))=- — In [(z - w)(z - w)] , 

(Uz, z)6j(w, vv)> = - J£ In [(z - w)(z - w)] , (2 ' 9) 
(fr(z,z)dj(w,w)} =0. 

Imaginary time ordering is implied for all correlation func- 
tions here and throughout the paper. Different bosonic fields 
above are of course uncorrelated. However, as we will show 
below, the presence of a junction mixes these fields and effec- 
tively reduces the independent bosonic degrees of freedom by 
half. 

At the end of wires, x = 0, fermions can hop between dif- 
ferent wires, this process is described by a single particle hop- 
ping Hamiltonian 

Hb = ~ Z M" </'.;<0RM» + H.c] , (2.10) 

U 



3 



where r, y and ary are the strength and the phase of the hop- 
ping amplitude between wires ; and j. Without loss of gen- 
eralit y 22 ' 27 (at least for junctions of two or three wires), we 
only consider the symmetric case fy = t in this paper. The 
phases encode the distribution of magnetic fluxes at the junc- 
tion, which for a junction of three or more wires can play a 
crucial role in the properties of some RG fixed points^ 2 ^ 



III. JUNCTION OF TWO WIRES: DEBC ANALYSIS 

In this section, we analyze the stability of RG fixed points, 
and compute their corresponding conductances for a junction 
of two TLL quantum wires with unequal Luttinger parame- 
ters. By using the DEBC method of Ref. |27l we show that 
the properties of such junctions only depend on an effective 
Luttinger parameter 

g; 1 = fa 1 + g?)n. (3.D 

Briefly, the junction of two wires is totally decoupled for 
g e < 1, and has a conductance g e (e 2 /h) for g e > 1. This 
result is consistent with what was found in Ref. [43[ where the 
tunneling between fractional quantum Hall edge states with 
different filling fractions was discussed. The DEBC method 
used in this paper, however, has the advantage that it can be 
straightforwardly generalized to junctions of more than two 
wires, (see SecfVland AppendixlAl 



A. DEBC method 

For one-dimensional quantum impurity problems, one of- 
ten invokes the conformal symmetry of the bulk system (Lut- 
tinger liquids in our case) and assumes that the effect of an 
impurity (junction), at low energies, is imposing a confor- 
mally invariant boundary condition (BC), which describes 
the renormalization-group (RG) fixed point. This method- 
ology of relating the BC and RG fixed points of the system 
is called boundary conformal field theory (BCFT), and has 
proved greatly successful in the study of quantum impurity 
problemsi^i 

A useful technique within the framework of BCFT is the 
DEBC method, which hugely simplifies the evaluation of the 
scaling dimensions, Aq b , of boundary operators, Ob, with a 
given BCiS The scaling dimension, in turn, determines the 
leading scaling behavior of a given operator under the RG 
flow, and thus govern the stability of the RG fixed points: In 
general, an RG fixed point (boundary condition) is stable if 
all boundary operators are either equivalent to identity or ir- 
relevant A<9 8 > 1 . Moreover, the conductance associated with 
the given fixed point can be readily computed from the BC. 
For a complete description of the DEBC method, we refer the 
reader to Refs. l27j|29l Here, we simply apply this method to a 
junction of two quantum wires with unequal Luttinger param- 
eters, commenting only on some key ingredients. 



B. Junctions of two quantum wires 

The first step of the DEBC method is to write an ansatz for 
the conformally invariant BCs describing the RG fixed points. 
The next step is to list all the boundary operators, which can 
possibly become relevant and make the fixed point unstable, 
and compute their scaling dimensions with such ansatz for ev- 
ery point in the parameter space, g\ and go- If for a given 
ansatz, none of these boundary operators have a scaling di- 
mension smaller than one (in some region of the parameter 
space known as the stability region), we have found a stable 
RG fixed point. In case of a junction, a natural ansatz can be 
expressed in terms of a rotation matrix ft that relates outgoing 
to incoming bosonic fields: 

s = <R$ L , (3.2) 

where cj> L ' R = (cf>^' R , . . . , (j^ ,R ) T are /-component vector fields. 

The most important boundary operators, in case of a junc- 
tion of two wires, correspond to the following processes: tun- 
neling of chiral fermions between the two wires, and backscat- 
tering within the individual wires. It is useful to introduce 
a compact notation for boundary operators describing the 
single-particle tunneling processes, 

T)° = <A -VfU, (3.3) 

where a,b — R,L. We can then list these six fundamental 
boundary operators in terms of the rescaled boson fields as 
follows: 

T n ( 22) -e'^VFS.m (3 4) 

j,ll{rr) i(J^-JjL) ±^(VFei-VH«a) 

21 e - e - 

The boundary operators corresponding to multiparticle pro- 
cesses are not forbidden and can be generated as higher-order 
perturbation processes even they are not presence in the bare 
Hamiltonian. In general, they can be constructed from these 
fundamental boundary operators and have larger scaling di- 
mensions and are less relevant than the single-particle pro- 
cesses. 

All the above boundary operators have the generic form 
Ob ~ e ,0 < * ' + ' be , where a, b are vectors that contain the pref- 
actors of the £>, and §, fields. By eliminating the redundant 
degrees of freedom with Eq. (13.2b . and using Eq. (12.9) , the 
scaling dimension of the generic Ob above can then be writ- 
ten in terms of % as 

= l -\<R T {a-b) + {a + b)\\ (3.5) 

where the superscript T represents matrix transpose. 

To find all the % matrices that correspond to stable fixed 
points, it is convenient to express £>, fields in terms of the fol- 
lowing ± fields: 

<P+ = ; = — ; <P- = ; = — • (3-6) 
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Corresponding 6 ± are defined in a similar manner. The six 
fundamental boundary operators can then be written as 



RL 

21(12) 
RL 

11(22) 



1/ i 



(3.7) 



r LL(««) <ig= iiVfefl- 
2 21 



Here, we have dropped all e' e * terms as they are effectively an 
identity at the boundary: Charge conservation requires 
j\ = 0, which, using Eqs. (12.51 and (12.61 l. gives <9 T #+ = (i.e., 
Dirichlet BC on 6+), and makes e ,fl+ an effective identity^ The 
simplified boundary operators in Eq. (13. 71 then only depend 
the effective Luttinger liquid parameter g,,. 

In terms of the left- and right-moving ± fields defined in a 
similar manner to Eq. (12.7b , the Dirichlet BC on 6+ gives 



A-0" 



Now we only need to specify the BC relating the iR fields 

as: 



Ia-=0 - 



(3.9) 



Because there is a single pair of RL fields, only the Neumann- 
BC (N-BC) and the Dirichlet-BC (D-BC), Kl' D = ±1, are 
allowed. By using Eq. ( 13. 5t on boundary operators listed in 
Eq. (13. 7t , we obtain the scaling dimension of each operator 
with the N-BC and D-BC: 



B 


A£ b (N-BC) 


A* (D-BC) 


T RL T RL 
1 21 ' 1 12 







T RL T RL 
1 1 1 ' 1 22 





& 


T LL T RR 

1 21 ' 1 21 


Vge 


ge 



Here, the scaling dimension indicates that Ob is equivalent 
to identity operator 1 for the given boundary condition. 

From the table above, we conclude that the N-BC is sta- 
ble when g e < 1 and D-BC is stable when g e > 1. As 
shown below, the N-BC corresponds to a fixed point where 
two wires are disconnected and the D-BC corresponds to a 
fixed point where two wires are maximally connected with 
conductance Go = ge(e 2 jK). Let us recall the well-known re- 
sults of the Kane and Fisher problem, namely, a single impu- 
rity in a spinless TLL wire, which is equivalent to a junction 
of two wires with the same Luttinger parameter, g. 14 There, 
two RG fixed points are identified, totally disconnected fixed 
point for g < 1 and maximally connected fixed point for g > 1 
with Go = g(e 2 /h). As g e — > g when g\ — g2 = g, the fixed 
points, which we identified, generalize Kane and Fisher's re- 
sults to the case of two quantum wires with different Luttinger 
parameters. 

A similar DEBC analysis can also be applied to junction of 
two quantum wires with spin-1 /2 electrons^ In that case, we 



obtain two effective Luttinger parameters g e 



2/(i/s r + 



l/g,' ), corresponding to charge and spin sectors, which gov- 
ern the scaling behavior of boundary operators near the fixed 
points^ As a consequence, the stable fixed points of such sys- 
tem are connected to those of equal Luttinger parameters but 
withg e "'^^' ! »i£ 



C. Conductance for each fixed point 

Here, we will compute the conductance associated with N- 
BC and D-BC. The conductance tensor is defined through the 
following current-voltage relation: 



(3.10) 



where current is defined as positive when flowing toward junc- 
tion. In the linear-response regime, the conductance above can 
be evaluated via the Kubo formula: 27 



dj = lim 

a)— >0+ 

(3.8) ~JZlJ_ dTe>b " X dx(Ji(y ' T)Jj(X ' 0)> ' (3 ' 1 1} 



The current correlation function can be rewritten as a sum of 
chiral-current correlation functions as 

<JiCy,Tyxjc,o)> = 

(Jf(y, r)jf(x, 0)> + (Jf(y, T)jfa, 0)> 
- (Jf(y, r)jf(x, 0)> - (Jf(y, t)J«(x, 0)), (3.12) 

where the left and right currents are defined in Eq. (12. 5t . By 
using Eq. ( 12.61 . we can write the chiral currents in terms of 
rescaled boson fields as follows: 



V2tt ~ 



V2tt 



(3.13) 



V27T V2tt 



where we have defined left and right currents associated with 
the rescaled fields. 

To evaluate the above correlation functions given a BC, it is 
convenient to first express the boundary conditions fi directly 
in the rescale boson field, cf>f' R , basis. From Eq. (13 . 8b and 
Eq. d3.9l l, we can derive the rotation matrix: 



2 *Jgi& 



g\ +gi 
gi-gi 

V gi+g2 8l+«2 



(3.14) 



for the N-BC and D-BC, respectively. From Eq. (13.131 . the left 
and right currents are constrained by the Jf (0, r) = flijJ^iO, r) 
BC at the origin, which, upon unfolding the current, implies 



J?(x,T)=%^(-Jt,T), 



(3.15) 



for x > 0. Therefore, all the right-moving currents can be 
interpreted as left-moving on the x < domain. 

Now, the chiral current correlation functions can be evalu- 



ated: 



(J?(Zi)J«(Zj)) - 



4tt 2 - zj) 2 ' 

(Jf(zdJj(Zj)) 8i . 2 , 

1 4n 2 {zi-zj) 2 



(3.16) 



(Jt(Zi)jf(Zj)) 



4n 2 (zi-zj) 2 ' 



By inserting these correlation functions into the Kubo formula 
Eq. (13. lit , and after some algebra, we obtain a concise rela- 
tion between conductances and boundary conditions^ 



(3.17) 



With the N-BC and D-BC represented by the rotation matrices 
of Eq. ( I3.14l i. we immediately conclude that 



G N = 0: 



1 -1 
-1 1 



(3.18) 



As expected, the N-BC corresponds to a fixed point with de- 
coupled wires and D-BC corresponds to a fixed point with 
conductance g e (e 2 /h). 

Here, we shall emphasize that the correlation functions 
listed in Eq. ( 13.161 1 include only the universal part for a given 
boundary condition. There are also nonuniversal contributions 
to the correlation functions, which, in general, decay faster 
and become irrelevant at long distance. However, when the 
universal part vanishes, which is the case for correlations be- 
tween different wires with the N-BC, the higher-order contri- 
butions dominate and could lead to nonlinear conductance. In 
the next section, we will perform the DMRG calculations and 
confirm that Eq. (13 . 1 6b indeed represents the universal part of 
chiral-current correlation functions. 



IV. JUNCTION OF TWO WIRES: DMRG CALCULATIONS 

In this section, we perform numerical computations of the 
conductance of a junction of two Luttinger liquids with dif- 
ferent Luttinger parameters g[ and g 2 . A microscopic lattice 
model with the Luttinger-liquid physics, which is suitable for 
numerical calculations, is the one-dimensional tight-binding 
model of interacting spinless electrons, 



c im Ci, m+ \ +H.c. + V,(«,> 



2 )(«!./»+ 1 



\\ (4.D 



for wire i, where the hopping amplitude is set to unity. A 
junction of two wires can be described by 



H = Hi + H 2 - tc[ c 2 ,o - tc\ ci, . 



(4.2) 



The parameters g, and v, of the Luttinger-liquid Hamilto- 
nian ( 12. U are related to the interaction strength V through the 
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FIG. 2: a) The junction of two semi-infinite Luttinger-liquid wires 
with interaction strength V[ and V 2 . Hopping amplitude is set to 
unity in the bulk of the the two wires, and is equal to / at the junction, 
b) The corresponding finite system used for DMRG calculations ac- 
cording to the method of Refs&S. 



Bethe ansatz (see, e.g., Refs. I66ll69h . We use the method of 
Refs. 66 67 to compute the conductance. This method allows 
us to extract the conductance from a ground-state static calcu- 
lation in a finite system as explained below. The semi-infinite 
junction of Hamiltonian ( 14.21 ) is depicted in Fig. IHa), while 
the corresponding finite system used in the numerics is shown 
inFig.Hb). 

We extract the conductance Gy from the following asymp- 
totic (large x) relationship in the finite system shown in 
Fig.Hb): 



ir) 



(4.3) 



where (J' R (x)J J L (x)) is the ground-state correlation function of 
chiral currents in wires i and j in Fig. [2t>- In terms of total 
charge and current / operators, which can be modeled on 
the lattice, we generically have 67 



(f R (x)J J L (x)) 



1 



2viV 



■{J\x)j\x)) ■ 



1 

27 



(N'(x)J J (x)). (4.4) 



Note that in the time-reversal symmetric case considered here, 
the second term in the expression above vanishes, and we only 
need to compute a static current-current correlation function 
{J'{x)jj{x)). In terms of the lattice creation and annihilation 
operators appearing in Hamiltonian (14. IK we have 



J'(m + -) = i(c 
2 



All we need to do now is to numerically compute (J 1 (x)J 2 {x)) 
for the above current operator, and divide it by 2 vi v 2 to obtain 
(7^(x)7^(x)). The numerical calculations are done for system 
of 180 sites in each of the two wires. The truncated number 
of states in our DMRG computations is 1 100. 



A. Repulsive effective interaction g e < 1 

For g e < 1, we have G\ 2 = 0, which implies that the 
leading term [Eq. (14. 3H in the (J^(x)J^(x)) correlation func- 



tion vanishes. If, as a function of | sin(|x), the computed 
(J^(x)J^(x)) decays faster than a power law with exponent -2, 
we have a signature of a vanishing G\o- Our numerical results 
indeed confirm this: for any combination of g\ and g2 with 

g e < 1, we find that {J x R (x)J^{x)) decays as ^sin^jc) ' 
with a(g e ) > 2. The exponent a only depends on g e , and is 
independent of individual g, and the hopping amplitude t . The 
prefactor of the correlation function depends on the hopping 
amplitude t since universality is a property of the leading term, 
and the coefficient of the subleading term observed here can 
depend on microscopic details such as t. 

Note that the correlation functions for different combina- 
tions of g\ and g2, which have the same g e , collapse not only 
in the large x limit but also close to the microscopic length 
scales. This strongly suggests that the single parameter g e 
determines all the subleading corrections to the correlation 
function. This behavior can be understood by noting that all 
the boundary operators in Eq. (13. 1\ depend on g e (as opposed 
to individual gi and g2) after dropping the tp + and 6 + fields 
due to current conservation. Hence, the measured correlation 
functions should also be determined only by the effective Lut- 
tinger parameter g e and the hopping strength t. Interestingly, 
the exponent a is indeed close to the scaling dimension of the 
leading irrelevant operator (i.e., 2/g e ) but we are not able to 
make a definitive statement due to finite-size effects and lim- 
ited numerical precision. It is worth mentioning that we have 
also considered combinations of g\ and g2 where one wire has 
attractive interactions {gi > 1), but this does not affect the be- 
havior of the junction as long as g e < 1 . 

The results are summarized in Fig. [3] We have considered 
two values of g e = 0.871,0.83 and four combinations of g\ 
and g2 for each g e as shown below. 
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Here, we have considered two different values of hopping am- 
plitude, t = 0.5,0.7. 
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FIG. 3: The correlation function {J l R (x)J 2 L (x)} for different combina- 
tions of gi and g 2 corresponding to two values of g e < 1. The data 
exhibits very good collapse, and indicates a vanishing conductance. 
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Again, we obtain very good agreement with the analytical pre- 
dictions. The results are summarized in Fig. |4] With fixed g e , 
the data points collapse for different combinations of g\ and 
g2 and different values of hopping t. The numerical results ap- 
proach the analytical prediction for a conductance of g e e 2 /h 
(solid black line) in the asymptotic limit. 



V. JUNCTION OF THREE WIRES 



B. Attractive effective interaction g e > 1 

In this case, we expect any tunneling amplitude t + 1 at the 
junction to heal, and result in a universal conductance gee 1 /h. 
In our numerics, we have used two values of tunneling ampli- 
tude, t = 0.7, 0.9. Similarly to the repulsive case, we consider 
two values of g e = 1.175, 1.258; for each g e , we consider four 
combinations of gi and g2 shown below. 



A junction of three quantum wires with equal Luttinger pa- 
rameters has three distinct types of fixed points described by a 
rotation matrix ansatz: decoupled fixed point, chiral-^ ± fixed 
points, and Dirichlet fixed points, which are respectively sta- 
ble for g<\, \ <g<3 and g > 3. (There is an addi- 
tional less understood time-reversal-invariant M fixed point 
for 1 < g < 3, which we do not consider in this work.) 22 ! 27 ' 66 ' 67 
Thus, it is reasonable to expect that these fixed points remain 
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FIG. 4: The correlation function (Jj i (x)Jl(x)) for different combina- 
tions of gi and g 2 corresponding to two values of g e > 1. The data 
exhibits very good collapse, and indicates a conductance of g e e 2 /h 
because of the asymptotic agreement with the exact theoretical pre- 
diction shown by a sold black line. 



stable if we slightly change the Luttinger parameters and make 
them unequal. Here, we first determine this region of sta- 
bility around the g\ = gi = gs line. In addition, we iden- 
tify three asymmetric fixed points, only realized for unequal 
Luttinger parameters, in which one of the wires is decoupled 
from the junction and the other two wires are fully connected. 
These asymmetric fixed points have important consequences 
for the stability of N, D, and^ ones: there are regions of the 
parameter space near the transition points gi = 1 , 3 on the 
gi - gi — g3 line where small perturbations normal to the 
equal-g line would drive the system into one of the asymmet- 
ric fixed points. 

In Table U we first summarize the scaling dimensions of 
the leading boundary operators as well as their corresponding 
conductances as a function of the three Luttinger parameters 
for each boundary condition. These scaling dimensions deter- 
mine the stability of the fixed point: when all them are larger 
than one in certain parameter region, the given fixed point is 
stable. From Table|I] we observe that the three scaling dimen- 
sions of the leading irrelevant operators for the asymmetric 
fixed points are each the inverse of those of the decoupled, 
Dirichlet and chiral fixed points. Hence, in any given point of 
the (gi , g2,g3) parameter space, there exists at least one stable 
fixed point. In other words, the decoupled, Dirichlet, chiral 
and asymmetric fixed points fully cover the phase diagram of 
a Y junction of spinless TLL wires. 

In this paper, the DMRG analysis is not applied to the junc- 
tion of three quantum wires with unequal Luttinger parame- 



ters. Such analysis with equal Luttinger parameters was per- 
formed in Ref. [63, but is beyond the scope of the present pa- 
per. We discuss the physical properties of each of these sta- 
ble fixed points in the remainder of this section. The detailed 
analysis, based on the DEBC method, can be found in Ap- 
pendix [A] As the conductances of fixed points can be evalu- 
ated in the similar way as in Sec. IIII CI we simply write down 
the results in the following discussion. 



a. Decoupled fixed point 

The decoupled fixed point corresponds to the Neumann BC 
for all the bosonic fields (p. From TableQ] we see that the scal- 
ing dimensions of the leading irrelevant operators are equal to 
1 Ig'e for g'J = 2gigj/(gi + gj), which is the same as Eq. d3.lt 
for a pair of wires. Hence, the decoupled fixed point (N-BC) 
is then stable when the N fixed point is stable for all three pos- 
sible pairs of wires. One can simply check that these scaling 
dimensions reduce to 1 /g when the Luttinger parameters are 
all equal. In Fig.(5j the stability region of the decoupled fixed 
point, A^J > 1, is shown in red. As expected, the conductance 
of the decoupled fixed points is simply 



G^ = 0. 



(5.1) 



b. Dirichlet fixed point 



As discussed in Appendix lAl one can construct three inde- 
pendent linear combinations of the bosonic field (p such that, 
akin to ip + field for junctions of two wires, one of them, known 
as the center of mass field, always satisfies the Neumann BC 
due to charge conservation. The Dirichlet fixed point corre- 
sponds to imposing the D-BC on the other two combinations. 
None of single particle processes becomes identity with such 
boundary condition. Instead, some of two- or more-particle 
processes reduce to identity, which suggests that the Dirich- 
let fixed point is associated with a certain type of "Andreev" 
reflection that enhances the conductance. 

From Table [I] the scaling dimensions of all leading irrele- 
vant operators, Ag s = g { (gj + g k )l2(g\ + g 2 + 83) V i * j + k, 
reduce to g/3 when gi = g2 = g3 — g- Hence, the D-BC 
becomes stable at g > 3, consistent with Ref. |27l By requir- 
ing > 1, we obtain the stability region, shown in green 
in Fig.ls] of the Dirichlet fixed point for unequal Luttinger 
parameters. 

The conductance of the Dirichlet fixed point is given by 



G * = 2 T 



gjgk 



g\+g2 + g3 



(5.2) 



where we have made a cyclic identification go = g3 and 
g 4 = gi. When all the Luttinger parameters are equal, we 
have G = g(e 2 /h)(26jk - 2/3), which reproduces the result 
in Ref. i3 

c. Chiral-^ ± fixed points 
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TABLE I: The scaling dimensions of leading irrelevant boundary operators and the conductance tensor for each stable fixed point of a Y 
junction. The detailed analysis for obtaining these scaling dimensions is given in Appendix lAl with the corresponding operators listed in 
Table Ull The asymmetric fixed point A, represents a boundary condition where the wire i, for ( = 1,2,3, is decoupled from the junction. 
Here, the scaling dimensions of all leading irrelevant operators run over the indices for all possible combinations of i + j ± k. We have also 
introduced following notations: the cyclic identification g = g 3 and g 4 = g, ; the two indices antisymmetric tensor e, ;±1 = ± 1 and otherwise; 
and the index m satisfies m + ) 4- k. The conductance tensors are given units of e 2 /h and are defined through Ij = Yik Gj k Vk- 



Fixed point 


Scaling dimensions A 


Conductance G s t [e 2 /h] 


Decoupled (N-BC) 


(gi + gj)/2gigj 





Dirichlet (D-BC) 


giigj + g*)/2(gi + gi + g 3 ) 


2[gj5jk - gjgkl(g\ + gl + g3>] 


Chiral (y±-BC) 


2 gi(gj + gk)l(glglg3 + g\ + gl + g3) 


T Sj (S 1 +S2 +S3 V>jk +gjgk( *gm ejk - 1 ) 
«1«2«3+.?1+S2+S3 


Asymmetric (A r BC) 


(gl + gl + g3 + glg2g3)/2gi(g i+l + gi-,) 

2gi+igi-i/(gi+i +gi-i) 

2(gl + gl + g?)lgi(gi+\ + gi-l) 


2fl|tfu.(-l + S tj + S it + 26 jt - 3<S„<5,,) 



The chiral-^ ± fixed points have a particular transport fea- 
ture: the realization of^- + or^f_ fixed points, with the incom- 
ing current respectively flowing clockwise or counterclock- 
wise into one of the adjacent wires, depends on the direction 
of the threaded magnetic field into the ringi^ (This point will 
become more apparent in the next section when discussing 
Fermi-liquid leads) When g\ — g2 — g3 = g, the scaling di- 
mensions listed in Table J] for both \± fixed points reduce to 
j^, and hence x± fixed point is stable for 1 < g < 3, which is 



consistent with what found in Ref. |27l Again, we obtain the 



stability region of the chiral fixed points through A*,* > 1 . In 



Fig. [5] such region is shown in orange. 

The conductances for chiral-^ ± fixed points are, in turn, 
give by 



G x ± _ 2 <? 2 gjigl +82 + 83)Sjk + gjgk(+gme jk - 1) 
jk h glg2g3 +gl +g2 +g3 

where e jJ±i = +1 while e jk - Ofor; = k, andg,„ + gj,g k . 



(5.3) 



d. Asymmetric fixed points 

The asymmetric-A,- fixed points have the property that wire 
i is decoupled from the junction. Such fixed points are a new 
feature of a system with unequal Luttinger parameters. As 
shown in Table Q] and mentioned earlier, the scaling dimen- 
sions of the leading irrelevant operators of A, fixed points are 
inverse to those of the decoupled, Dirichlet and chiral fixed 
points. Hence, the stability regions of A, fixed points are com- 
plimentary to those of the other fixed points, and the entire 
parameter space is covered by at least one of the stable fixed 
point presented in this work. In Fig. [5] the stability regions of 
Ai 2,3 fixed point are shown in yellow, gray, and blue, respec- 
tively. The regions where two asymmetric fixed points overlap 
are shown in white. 

The conductances of the asymmetric fixed points A, are give 

by 



G A i = — (-1 + Sij + S a + 26 jk - 35y5 ft ), 



'jk 



(5.4) 



where g"'' n is the effective Luttinger parameter for the pair of 
wires m and n. To give an idea of the properties of asymmetric 



fixed points, we hereby write the conductance tensor explicitly 
for Ai, 



2g2g3 e 2 



g2 + g3 h 



(0 
1 
-1 







1 



(5.5) 



The stability of the different fixed points in the g 1,2,3 pa- 
rameter space is inferred from the scaling dimensions of Ta- 
ble U and shown in Fig. [5] Since each stable fixed point im- 
plies a phase of the junction, we shall also refer to the graph 
defining the regions of stability as a phase diagram. It is con- 
venient to illustrate this three-dimensional phase diagram by 
some cross sections. As all the scaling dimensions of lead- 
ing irrelevant operators have a cyclic symmetry on wire in- 
dices, the stability regions show a threefold rotation symmetry 
around the g\ - g2 = g3 axis. Thus, a natural choice for these 
cross sections is given by planes normal to the equal-g axis: 
(gi +g2 +g3)/3 = g, where the parameter g, the average of the 
three Luttinger parameters, labels each cross section. Exclud- 
ing negative gi, these cross sections are equilateral triangles 
shown in Fig. [5] 

We first notice that the decoupled fixed point becomes pre- 
dominant when g < 2/3, in which none of other fixed points 
are stable. One can also show that chiral fixed points appear 
when g > 1 and the Dirichlet fixed point only appears when 
g > 3. From Figs. [3b), [2d), and|3f), we observe that the de- 
coupled, chiral and Dirichlet fixed points are realized around 
the equal-g axis for g < 1,1 < g < 3 and g > 3, respectively. 
These results are consistent with a junction of three identical 
TLL wires (indicated as black points at the center of triangles), 
and show how far the Luttinger parameters can deviate from 
equal-g axis before these phases break down. We find that, 
in most regions, these phases are stable under a small pertur- 
bation away from the equal-g line. This is relevant for the 
realization of these phases experimentally, as the TLL wires 
attached to a junction are likely nonidentical. 

As shown in Figs.|5jc) and|5je), the asymmetric fixed points 
become important around two marginal points £1,2,3 = 1,3. 
As a small deviation of Luttinger parameters from these two 
points easily realizes and switches between asymmetric fixed 
points, the resultant fixed points are thus highly sensitive to all 
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FIG. 6: The voltages for the incoming and outgoing chiral currents 
in the Fermi-liquid and connected TLL wire. 



the three Luttinger parameters, and not just the averaged Lut- 
tinger parameter g. Therefore, precise control over the TLL 
wires become essential around these points. Note that it may 
be possible to alter between different asymmetric fixed points, 
and form a nano-switch if one can tune the Luttinger param- 
eters. We mention in passing that the Luttinger parameters of 
wires can, in principle, be modified by an external gate capac- 
itively coupled to the wire4i 

Finally, it is worthwhile to compare our findings with those 
of Aristov and Wolfle^Si where two identical wires are con- 
nected to a wire with unequal Luttinger parameter, by set- 
ting gi = g2 — g- In the repulsive and weak attractive in- 
teraction regime, g m g3 < 3, in which the N, x± an d A,- 
fixed points predominate, our results show excellent agree- 
ment with their findings. In the strong attractive interaction 
regime, g « g^ > 3, the D fixed point identified in the present 
work and in Ref . |27] (by nonperturbative boundary-conformal - 
field-theory methods), however, was not found by the pertur- 
bative renormalization-group approach of Refs. P8ll59l 



FIG. 5: We plot the stability regions of the different fixed points 
in the triangular cross sections, shown in upper panel, of fixed 
g = (gi+^2+ft)/3 in Luttinger parameters space, g lt g 2 , and #3. In all 
panels, the decoupled, chiral and Dirichlet fixed points are painted, 
respectively, in red, orange, and green, and, the asymmetric fixed 
points A [ 2.3 are pained in yellow, gray, and blue, respectively. The 
white areas represent an overlap of two Asymmetric fixed points and 
the equal-g points are indicated by a black dot at the center of trian- 
gle, (a) For g < 2/3, only the decoupled fixed point is realized, (b) 
The^± and Dirichlet fixed points are not stable when g < 1. (c) The 
point gi 1,3 = 1 is an exactly marginal point surrounded by asymmet- 
ric fixed points, (d) The^ fixed points appear at the center of cross 
section for 1 < g < 3, surrounded by asymmetric fixed points. The 
stability region for the N-BC are pushed to the corners where two of 
the three Luttinger parameters become much less than one. (e) The 
point gi 2,3 = 3 is another exactly marginal point, again surrounded 
by asymmetric fixed points. The chiral fixed points have extremely 
small stability regions (difficult to see in the figure) located between 
any two asymmetric fixed points, (f) The Dirichlet fixed point ap- 
pears at the center of triangle for g > 3. Note that the figure is not a 
schematic, and represents the exact domain boundaries for the shown 



VI. CONDUCTANCE RENORMALIZATION FOR WIRES 
CONTACTED TO FERMI-LIQUID LEADS 

The linear conductances of different fixed points were cal- 
culated in Sees. IIII CI and [V] Here, we discuss the effect of 
attaching the wires to Fermi-liquid leads. Remarkably, we 
find that the conductance of each fixed point, in the presence 
of Fermi-liquid leads, renormalizes to values that are indepen- 
dent of the Luttinger parameters. This generalizes the follow- 
ing interesting effect for the TLL quantum wire with Luttinger 
parameter g: When attached to leads, the measured conduc- 
tance is quantized at e 2 /h, which is different from ge 2 /h. This 
discrepancy has been resolved by Maslov and Stone and by 
Safi and Schulz in Refs. 38 39J There, they studied an inho- 
mogeneous Luttinger liquid and concluded that the conduc- 
tance of a TLL wire will only depend on the Luttinger pa- 
rameter at the contact. As a Fermi-liquid (metal) contact can 
be thought of as a TLL with Luttinger parameter g = 1 , the 
measured conductance becomes simply e 2 /h. 

An alternative way to understand this renormalization of 
conductance due to a Fermi-liquid contact is to introduce a 
contact resistance, 1/G C = (g - l)(2g)(h/e 2 ), at both contacts 
in series with the theoretical predicted resistance l/g(h/e 2 ). 
This would give the total conductance exactly at e 2 /h. Here, 
the contact conductance causes a voltage drop when matching 
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with the wire and gives rise to a current/voltage relation 

k = G\{Vi - Vj), (6.1) 

where V, is the potential of the carriers injected into wire i and 
Vj is the applied voltage at the contact connected to the wire 
i. We can understand this result by thinking in terms of the 
right- and left-moving currents inside the wire and the Femi- 
liquid lead as seen in Fig. [6] The incoming electrons from 
the Fermi-liquid side are at the voltage V, of the reservoir. 
The outgoing current may be at a different potential V° at at 
the contact point, even though the electrons are expected to 
relax to the equilibrium voltage as they propagate in the lead. 
Similarly, the right- and left-moving currents in the TLL wire 
are respectively at voltage V, and V™. The junction of the 
Fermi-liquid lead and the TLL has an effective conductance 
g e e 2 /h with 

l/& = (l + l/*)/2. 

The current is then related to the difference in the voltages of 
the two incoming currents as 

h = geiVi - V] n ). 

For the TLL wire, we can similarly write 

ii = g(Y, - v'T). 

Combining the above two equations leads to Eq. ( 16.11 ). Gen- 
erally, it is useful to define a contact conductance tensor 
G c = 5ij(e 2 /h)(2gj)/(gi - 1) for a junction of three wires. 

As the measured conductance is based on the applied volt- 
age at contact, one can define a renormalized conductance ten- 
sor, Gij, as I, = Yij &ijVj. To connect this conductance with 
one we found in the previous subsection, we invoke current 
conservation: 

4 = 2>; y ; = Z^- (6 - 2) 

j i 
Together with Eq. ( I6.lt . one can show that^ 

G = {\ + GG- l y l G, (6.3) 

or equivalently 

G 1 = G 1 + G~\ (6.4) 

which has a simple interpretation of resistances connected in 
series. 

Besides the decoupled fixed point that has obvious vanish- 
ing conductance G N = 0, we shall now apply Eq. d6.3t and 
obtain the measured conductances of the other fixed points. 
For Dirichlet fixed point, we have the renormalized conduc- 
tance 

<3g = ~(35y-l). (6.5) 

As the largest conductance one can obtain from single-particle 
unitary scattering is G u , = —(4/9)(e 2 /h) for j + the 



enhanced conductance above demonstrates the role of multi- 
particle scattering processes. 

Upon attaching the wire to external Fermi-liquid leads, the 
measured conductances of chiral fixed points become 

^ = ^x[ (3 ^ _1)T ^] (6 ' 6) 

(i.e., the currents flow only from lead 1 to 2, 2 to 3, and 3 
to 1 for the x + fixed point and in reversed order for x fixed 
point). Finally, the measured conductance of asymmetric-A, 
fixed points reads 

-a e 2 

Gfj = j(-l + Stj + So, + 26 jk - 35ijS ik ), (6.7) 

which simply indicates a decoupled wire i with the rest of two 
wires fully conducting. 

We note that all renormalized conductance G D , &*, and 
G Ai , are the same as the unrenormalized conductance G with 
all gj = 1. This result highlights that the dc conductance 
of a junction of TLL wires depends only on the asymptotic 
value of the Luttinger parameters of the wires, and in the case 
when Fermi liquid leads are attached, this asymptotic value 
is gi = 1. Thus, when in contact to Fermi liquid leads, the 
conductance tensor of the junction will take on the universal 
values listed above for the different fixed points. Notice, how- 
ever, that which fixed point is selected still depends on the gj's 
of the interacting wire segments. 

VII. CONCLUSION 

In summary, we applied the DEBC method to a two-wire 
and a Y junction of TLL wires with generally unequal Lut- 
tinger parameters. For two spinless wires, we successfully re- 
produced the prediction that all properties of the junction are 
determined by a single effective Luttinger parameter g e . We 
verified this prediction by direct numerical calculations on the 
lattice with the method of Refs. l66ll67l : We observed numeri- 
cally that as long as g e < 1, even if one of the wires has at- 
tractive interactions gj > 1, any impurity leads to a vanishing 
linear conductance. Moreover, we found that the nonuniver- 
sal corrections to the correlations across the junction, which 
come from perturbations with irrelevant boundary operators 
to the decoupled fixed point, are independent of the individ- 
ual Luttinger parameters and only depend on g e and the local 
microscopic structure of the junction. For g e > 1, we explic- 
itly found a universal conductance of g e e 2 /h regardless of the 
individual Luttinger parameters and the microscopic details. 

For a Y junction of nonidentical TTLs, we found that the 
N-BC, ^-BC and the D-BC are stable within regions of the 
(g1.g2.g3) parameter space, which we explicitly determined. 
By identifying three more asymmetric fixed points, corre- 
sponding to only one decoupled wire, and determining the re- 
gion of stability of each fixed point, we determined the full 
phase diagram of the Y junction. We also obtained explicit 
formulas for the conductance of all fixed points. The findings 
of this work have direct experimental relevance. In particular, 
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our results shed light on the issue of connecting interacting 
TLLs to Fermi-liquid leads for the measurement of transport 
properties. Our work also provides an important theoretical 
extension of the well-known results on transport through junc- 
tions of two and three identical TLL wires. 

As an outlook, we finally discuss how to generalize the cur- 
rent method to a junction of N > 3 wires. The key to such 
generalization is to identify the possible fixed points via the 
rotation matrix fi. These rotation matrices are constrained 
by charge conservation and can, in general, be expressed as 
S 0(N - 1) matrices after eliminating the total-charge (center- 
of-mass) mode (cf., Appendix lAl. Then, one has to identify 
the corresponding N-BC, D-BC, chiral-like BC, and asym- 
metric BC and analyze their stability. One useful trick for 
identifying the possible stable fixed points is to utilize the fact 
that a stable BC would make certain boundary operators ef- 
fectively equal to identity, cf. the discussion in Appendix lAfc. 
Of course, with the number of wires increasing, the number of 
possible fixed points also increases and the analysis becomes 
more complicated. 



TABLE II: The boundary operators corresponding to single-particle 
processes at the Y junction. 

(a) ± cycle 
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Appendix A: Junction of three wires with unequal Luttinger 
parameters: DEBC analysis 

In this Appendix, we present the details of the DEBC anal- 
ysis of the stability of the following fixed points for a junc- 
tion of three quantum wires: decoupled, chiral-^- ± , Dirich- 
let, and Asymmetric-A, fixed points. The system is described 
by the action ( 12. lb with i = 1,2,3 and the hopping Hamil- 
tonian (12.10b with ary = y/3, where y is the magnetic flux 
through the ring at the junction. 

To simplify the notation, we drop the overhead tilde (~) 
symbol for the rescaled fields throughout this appendix. (All 
the bosonic fields are rescaled.) To employ the DEBC method, 
it is convenient to choose a proper basis for the rescaled bo- 
son fields. In the first step, we identify the following center- 
of-mass field, which always satisfies the N-BC due to charge 
conservation: 



1 







y/gi +gi + g3 



The dual field to <l>o, i.e.. 

1 



Oo = 



Vli + 82 + g3 



( + V#2^2 + V§3^3)- (Al) 



( VSi&i + V^fe + VgJft) (A2) 



then becomes a constant and can be simply neglected. We 
then define another two orthonormal boson fields: 



3> 2 = 



y/gi + 82 
t y/gigm + y/g2gi<P2 



(gi +^2)^3] 



(A3) 



VSl +g2+g3 y/gi + g2 



as well as their dual fields, ©1,2- Note that the choice of ba- 
sis above is arbitrary but, as will become apparent, is a con- 
venient one. We will organize the fields above into a vector 
<j> = (d>i, 2 ) r and its dual vector = (©1, & 2 ) T . 
It is useful to define the following vectors: 



K x =(- 



K 2 =(~ 



igl + g2+g3 



y/2g2(gi + g2) y/2g 3 (gi + g 2 ) 

y[g~2 igl + g2+g3 



V 2 gi(gi + £2) ^2g 3 (gl + g2) 



). 



(A4) 



yj^gigi 

which will further simplify the notation. Notice that these 
three K , vectors add up to for any gj. 

By using the notation in Eq. ( 13.3b and neglecting ©o (due 
to the N-BC on <£>o), the boundary operators for single-particle 
processes are categorized in four classes, and listed in TablelHl 
As in the two-wire case, the higher-order processes can be 
constructed from these single-particle boundary operators. 



Using the ansatz (13.2b . we then write a rotation matrix: 



Re 



cos £ sin £ 
-sin£ cosf 



(A5) 



12 



where £ is a rotation angle. The rotation matrix above relates 
4> L = (* + 0)/2 to 0« = ($ - 0)/2. In terms of % the 
scaling dimension of boundary operators is given by Eq. (13.51 ). 
We now proceed to the stability analysis of the four Y-junction 
fixed points. 

a. Decoupled fixed point 

The decoupled fixed point corresponds to the N-BC for <& 
field, and makes the field a pure number. Therefore, all 
backscattering processes are effectively identity. The rotation 
matrix is simply equal to = 1 . From Eq. (13.5b , the scal- 
ing dimension of an arbitrary operator with the N-BC becomes 
A^ = \a\ 2 . The explicit scaling dimensions for the operators 
in Table im are listed below. 



E 



A% b (N-BC) 



T RL 
1 21 ' 


T RL 
1 12 ' 


T LL 
1 21 ' 


T RR 
1 21 


\K,\ 2 


= (gi 


+ 82)/(2gigi) 


T RL 
1 32 ' 


T RL 
1 23 ' 


T LL 
1 32 ' 


T RR 
1 32 


\K X \ 2 


= (gl 


+ gj)K2g2g-h) 


T RL 
1 13 ' 


T RL 
1 31 ' 


T LL 
1 13 ' 


T RR 
1 13 


\K 2 \ 2 


= (#3 


+ gi)/(2g3gi) 



Note that condition > 1 determines the stability region of 
the N-BC, shown in red in Fig. [5] 

b. Dirichlet fixed point 



where £ ± are the rotation angles of ^- ± -BC, respectively. The 
rotation matrices are obtained by plugging in the respective 
rotation angles into Eq. ( IA51 . Here, we do not show them 
explicitly. 

Let us first focus on the x+ fixed points. The scaling di- 
mensions of the single particle processes (excluding + cycle) 
are: 



B 



T RL 
1 W 
T RL 
1 22 ' 
T RL 
1 33 ' 



T RL 
1 23 ' 


T LL 
1 13 ' 


T 


T RL 
1 31 ' 


T LL 
1 21 ' 


T 


T RL 
1 12 ' 


T LL 
1 31 ' 


T 



RR 
21 
RR 
32 
RR 
13 



2gl(g2+gi) 
glg2g3+gl+g2+g3 

2g2<gl+g3) 
glg2g3+gl+g2+g3 

2g3(gl+gl) 
gjg2gl±Sl±g2+gj 



Notice that these scaling dimensions are cyclic in three indices 
and hence all operators are important for determining the sta- 
bility of the x+ fixed point. 

As for the^_ fixed point, one can show that all the leading- 
order operators have exactly the same scaling dimensions 
listed in the table above. Thus, both x± fixed points share 
exactly the same stability. In Fig. |5j the stability region of 
chiral-^ ± fixed point is shown in orange. 



d. Asymmetric fixed points 



The Dirichlet fixed point corresponds the D-BC on the 
<fr field (i.e., $ is effectively a constant at boundary). The 
rotation matrix of D-BC simply reads = -1. From 

Eq. ( 13.51 ), the scaling dimension of an arbitrary operator with 
D-BC becomes A^ = \b\ 2 . Unlike the N-BC, none of the 
single-particle processes becomes identity with D-BC. How- 
ever, some of the two-or-more-particle processes do become 



identity under D-BC, for instance: T^Tf^ , 



T RL T RL1 
1 32 1 23 



and 



Tf^T?^ . This indicates that the Dirichlet fixed point is asso- 
ciated with the Andreev reflection. 

Hereby, we only list the scaling dimensions for the ± cycle 
as they are the leading irrelevant operators. 



B 






A° b (D-BC) 


T RL 
1 21 ' 
T RL 
1 32 ' 
T RL 
1 13 ' 


T RL 
1 12 
T RL 

23 
T ~RL 
Ml 


gi{g\ 
g\{gl 
gz(&3, 


+ g2)/2(gl +g2+gi) 
+ g3)/2(gl +gl+gi) 
+ gl)/2(?l +g2+gi) 



Now, the Dirichlet fixed point is stable only when all these 
scaling dimensions A^ > 1. In Fig.[5j the stability region of 
Dirichlet fixed point (D-BC) is painted in green. 

c. Chiral-^ ± fixed points 

The chiral-^ ± fixed points are defined as follows: the BC 
corresponding to the x + (X ) fixed point would effectively 
make all operators in + (-) cycle equal to identity for all Lut- 
tinger parameters. With this in mind, we could derive the fol- 
lowing relationship for the rotation angles of the correspond- 
ing rotation matrices: 



tan£ ± = + 



gl + gl + g3 

glglg?, 



(A6) 



Although the use of rotation matrices fi^ is useful for other 
fixed points, it is most convenient to identify rotation matrices 
directly in the rescaled boson field, <^" s , basis. Because the 
decoupled wire effectively has the N-BC for itself and the con- 
nected wires should follow mutual D-BC, by using Eq. <I3 . 14b . 
the rotation matrix of asymmetric fixed point, A \ , has the form 



& = 



1 1 











g2~g3 


2 Vg2g3 


g2+g3 


g2+g3 


,0 


2 Vg2g3 


g3~g2 


g2+g3 


g2+g3 



(A7) 



while those of A23 can be constructed by permuting the in- 
dices in the corresponding matrix elements. 

By using this rotation matrix (IA7b . it is straightforward to 
show that the following single-particle tunneling processes are 
equal to identity: T^, T^t, and T^. In addition, the scaling 
dimensions of the leading relevant/irrelevant operators read as 
follows. 



B 




A* (Ax-BQ 


T RL T RL 
21 ' 12 ' 
T RL -rRL 
1 22 ' 1 33 ' 
T RL T RL^ 
1 21 1 12 ' 


T LL T RR T RL T RL T LL T RR 
1 21 ' 1 21 ' 1 13 ' 1 31 ' 1 13 ' 1 13 
T LL T RR 
32 ' 32 
T RL T RL^ 
1 13 1 31 


gl+g2+g3+glg2g3 
2gl(g2+g3) 
2g2g3 
g2+g3 

2(- + -M 

V fil g2+g3 7 



Here, we notice that some leading order operators are two par- 
ticle processes. To obtain the scaling behaviors of operators 
near the A23 fixed points, one can simply permute the indices 
of the Luttinger parameters with the corresponding operators. 
In Fig.|5j the stability regions of Aj 2,3 fixed points are shown 
in yellow, gray, and blue, respectively. 
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